Trajectory planning

ABSTRACT

A method and apparatus for determining a trajectory for a vehicle are disclosed, wherein, the method includes: identifying a starting position (p 0 ) and a desired terminal position (P) for the vehicle; linearly approximating dynamics of the vehicle; and using the starting position (p 0 ), desired terminal position (P), and linear approximation, determining the trajectory for the vehicle. The linear approximation can be constrained by requirements (e.g., specifications) that: (i) an acceleration applied to the vehicle at a point on the trajectory is relatively large when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; and (ii) an acceleration applied to the vehicle at a point on the trajectory is relatively small when the acceleration acts in a direction that is substantially parallel to the velocity of the vehicle. The vehicle may have a curvature limit.

FIELD OF THE INVENTION

The present invention relates to trajectory planning and apparatus for the planning of trajectories for vehicles.

BACKGROUND

Trajectory planning methodologies, for example using Mixed Integer Linear Programming (MILP), are used to determine globally optimal trajectories for vehicles. Many trajectory planning methodologies constrain vehicle trajectories with a linear approximation of the vehicle's dynamics.

Typically, a linear approximation of the vehicle's dynamics does not contain a notion of vehicle heading. This is typically because the introduction of a heading angle introduces non-linearities. As a result, a determined trajectory may feature motion that cannot be achieved by conventional wheeled or tracked vehicles.

SUMMARY OF THE INVENTION

In a first aspect, the present invention provides a method for determining a trajectory for a vehicle, the method comprising: identifying a starting position for the vehicle; identifying a desired terminal position for the vehicle; linearly approximating dynamics of the vehicle; and using the starting position, the desired terminal position, and the linear approximation, determining the trajectory for the vehicle; wherein the linear approximation is constrained by requirements that: (i) an acceleration applied to the vehicle at a point on the trajectory is relatively large when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; and (ii) an acceleration applied to the vehicle at a point on the trajectory is relatively small when the acceleration acts in a direction that is substantially parallel to the velocity of the vehicle.

The constraint that the acceleration is relatively large when acting in a direction that is substantially perpendicular to the velocity of the vehicle, and is relatively small when acting in a direction that is substantially parallel to the velocity of the vehicle, may be implemented using the following:

  ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circs)] ${{\left( {\frac{{a_{x}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{x}(k)}}{V_{{ma}\; x}}} \right){\cos\left( \frac{2\pi\; m}{N_{circs}} \right)}} + {\left( {\frac{{a_{y}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{y}(k)}}{V_{{ma}\; x}}} \right){\sin\left( \frac{2\pi\; m}{N_{circs}} \right)}}} \leq \sqrt{1 + \lambda^{2}}$ where: v_(x)(k) is the velocity of the vehicle in an x-direction;

-   -   v_(y)(k) is the velocity of the vehicle in a y-direction;     -   V_(max) is a maximum velocity of the vehicle;     -   N_(steps) is a number of time-steps for the trajectory;     -   N_(circs) is a number of corners of a polygon;     -   R_(min) is a minimum turn radius of the vehicle; and     -   λ is a constant.

The linear approximation of the vehicle dynamics may be further constrained by the requirement that a magnitude of a velocity of the vehicle is greater than or equal to a threshold value for the velocity. This constraint may be implemented using the following:

∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circ)] ${{{v_{x}(k)}{\sin\left( \frac{2\pi\; m}{N_{circ}} \right)}} + {{v_{y}(k)}{\cos\left( \frac{2\pi\; m}{N_{circ}} \right)}}} \geq {V_{m\; i\; n} - {{p\left( {k,m} \right)} \cdot C_{1}}}$ and ∀k ∈ [0, …  , N_(steps)] ${\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{circ} - 1}$ where: v_(x)(k) is the velocity of the vehicle in an x-direction;

-   -   v_(y)(k) is the velocity of the vehicle in a y-direction;     -   V_(min) is the threshold velocity value;     -   N_(steps) is a number of time-steps for the trajectory;     -   N_(circ) is a number of corners of a polygon;     -   p(k,m) is a binary decision variable; and     -   C₁ is a constant.

C₁ may be relatively large compared to V_(min).

The linear approximation may be further constrained by a requirement that the vehicle may not travel in a given region; a centre point of the region is at a distance from a predetermined point on the vehicle substantially equal to a minimum turn radius of the vehicle, in a direction substantially perpendicular to a velocity of the vehicle; and a distance from the centre point of the region to a point on a perimeter of the region is greater than or substantially equal to the minimum turn radius of the vehicle.

The constraint that the vehicle may not travel in a given region may be implemented using the following:

  ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(regions)] ${{\left( {{r_{x}(k)} - r_{0x} - {R_{m\; i\; n}v_{nx}}} \right){\cos\left( \frac{2\pi\; m}{N_{regions}} \right)}} + {\left( {{r_{y}(k)} - r_{0y} - {R_{m\; i\; n}v_{ny}}} \right){\sin\left( \frac{2\pi\; m}{N_{regions}} \right)}}} \geq {R_{m\; i\; n} - {C_{2}{p\left( {k,m} \right)}}}$   with   ∀k ∈ [0, …  , N_(steps)] $\mspace{20mu}{{\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{regions} - 1}}$ where: v_(nx) is the normalised velocity of the vehicle in an x-direction;

-   -   v_(ny) is the normalised velocity of the vehicle in a         y-direction;     -   r_(x)(k) is a position of the vehicle from an origin in the         x-direction;     -   r_(y)(k) is a position of the vehicle from an origin in the         y-direction;     -   r_(0x) is an initial position of the vehicle from an origin in         the x-direction;     -   r_(0y) is an initial position of the vehicle from an origin in         the y-direction;     -   N_(steps) is a number of time-steps for the trajectory;     -   N_(regions) is a number of corners of a polygon;     -   R_(min) is a minimum turn radius of the vehicle;     -   p(k,m) is a binary decision variable; and     -   C₂ is a constant.

The linear approximation may be further constrained by requirements that: an acceleration of the vehicle during the trajectory is less than a threshold value for the acceleration; and the threshold value for the acceleration is dependent on an infinity norm of a velocity of the vehicle. This constraint may be implemented using the following:

$\begin{matrix} {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{m \in \left\lbrack {1,\ldots\mspace{14mu},N_{cir}} \right\rbrack}}}{{{{a_{x}(k)}{\sin\left( \frac{2\pi\; m}{N_{cir}} \right)}} + {{a_{y}(k)}{\cos\left( \frac{2\pi\; m}{N_{cir}} \right)}}} \leq {a_{m\;{ax}}(k)}}} & (i) \\ {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{p \in \left\lbrack {1,\ldots\mspace{14mu},N_{pol}} \right\rbrack}}}{{{{m(p)}{v_{x}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 1\left( {k,p} \right)}}}{{{{m(p)}{v_{y}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 2\left( {k,p} \right)}}}} & ({ii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{\sum\limits_{p = 1}^{N_{pol}}\left( {{b\; 1\left( {k,p} \right)} + {b\; 2\left( {k,p} \right)}} \right)} \leq {{2N_{pol}} - 1}}} & ({iii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{a_{{ma}\; x}(k)} \geq 0}} & ({iv}) \end{matrix}$ where: a_(max) is a maximum acceleration of the vehicle;

-   -   a_(x) is an acceleration of the vehicle in an x-direction;     -   a_(y) is an acceleration of the vehicle in a y-direction;     -   v_(x) is a velocity of the vehicle in an x-direction;     -   v_(y) is a velocity of the vehicle in a y-direction;     -   N_(steps) is a number of time-steps for the trajectory;     -   N_(cir) is a number of corners of a polygon;     -   N_(pol) is a number of linear functions, each linear function         relating the maximum acceleration of the vehicle to an infinity         norm of the velocity of the vehicle;     -   m(p) is a value of a gradient of a pth linear function;     -   c(p) is a value of pth linear function if the value of the         velocity of the vehicle is zero;     -   b1(k,p) and b2(k,p) are constants; and

C_(L) is a constant.

The method may further comprise: determining one or more further linear approximations of the dynamics of the vehicle; and using the one or more further linear approximations, determining one or more further trajectories for the vehicle; wherein the trajectory and the one or more further trajectories are for implementation by the vehicle in series; each of the one or more further linear approximations is constrained by one or more of the following: (i) a requirement that a magnitude of a velocity of the vehicle is greater than or equal to a threshold value for the velocity; (ii) requirements that: an acceleration applied to the vehicle at a point on the trajectory is relatively large when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; and an acceleration applied to the vehicle at a point on the trajectory is relatively small when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; (iii) a requirement that the vehicle may not travel in a given region; a centre point of the region is at a distance from a predetermined point on the vehicle substantially equal to a minimum turn radius of the vehicle in a direction substantially perpendicular to a velocity of the vehicle; and a distance from the centre point of the region to a point on a perimeter of the region is greater than or substantially equal to the minimum turn radius of the vehicle; and (iv) requirements that an acceleration of the vehicle during the trajectory is less than a threshold value for the acceleration, and the threshold value for the acceleration is dependent on an infinity norm of a velocity of the vehicle.

The vehicle may have a curvature limit for its trajectory.

The vehicle may be a land-based vehicle.

In a further aspect the present invention provides apparatus for determining a trajectory for a vehicle, the apparatus comprising one or more processors arranged to: linearly approximate dynamics of the vehicle; and using an identified starting position for the vehicle, an identified desired terminal position for the vehicle, and the linear approximation, determine a trajectory for the vehicle; wherein the linear approximation is constrained by requirements that: (i) an acceleration applied to the vehicle at a point on the trajectory is relatively large when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; and (ii) an acceleration applied to the vehicle at a point on the trajectory is relatively small when the acceleration acts in a direction that is substantially parallel to the velocity of the vehicle.

In a further aspect the present invention provides a program or plurality of programs arranged such that when executed by a computer system or one or more processors it/they cause the computer system or the one or more processors to operate in accordance with the method of any of the above aspects of the invention.

In a further aspect the present invention provides a machine readable storage medium storing a program or at least one of the programs according to the above aspect of the invention.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a schematic illustration (not to scale) of a land based vehicle for which an embodiment of a method of performing path planning is to be implemented;

FIG. 2 is a schematic illustration (not to scale) of a scenario in which the vehicle travels along a trajectory;

FIG. 3 is a schematic illustration (not to scale) representing a velocity space of the vehicle;

FIG. 4 is a process flow chart showing certain steps of a process of determining a constrained acceleration vector for the vehicle;

FIG. 5 is a schematic illustration (not to scale) that is useful in the understanding of the process of FIG. 4;

FIG. 6 is a schematic illustration (not to scale) showing the vehicle and two regions into which the vehicle may not directly move; and

FIG. 7 is a schematic illustration (not to scale) of a graph showing the non-linear relationship between a maximum acceleration of the vehicle and an ∞-norm of the velocity of the vehicle.

DETAILED DESCRIPTION

FIG. 1 is a schematic illustration (not to scale) of a land based vehicle, hereinafter referred to as “the vehicle 2”, for which an embodiment of a trajectory planner (not shown) plans a trajectory. The terminology “trajectory” is used herein to refer to a route of the vehicle from a first position to a second position.

The following information about a state and operation of the vehicle 2, described with reference to FIG. 1 and FIG. 2, will be useful in understanding the embodiment described later below.

The vehicle 2 has the following state vector x:

$x = \begin{pmatrix} r_{x} \\ r_{y} \\ \theta \end{pmatrix}$ where: r_(x) is a distance of the vehicle 2 from an origin ◯ in the x-direction, as indicated in FIG. 1;

-   -   r_(y) is a distance of the vehicle 2 from the origin ◯ in the         y-direction, as indicated in FIG. 1; and     -   θ is an angle between the x-axis and a heading of the vehicle,         hereinafter referred to as “the heading” and indicated in FIG. 1         by a solid arrow and the reference numeral 4.

Also,

$\overset{.}{x} = \begin{pmatrix} {v\;\sin\;\theta} \\ {v\;\cos\;\theta} \\ \omega \end{pmatrix}$ where: {dot over (x)} is the time derivative of the state vector x;

-   -   v is a value of the speed of the vehicle that the vehicle         travels with (in the direction of the vehicle's heading 4); and     -   ω is the curvature of the vehicle's path.

In this embodiment, the vehicle 2 has the following input vector u:

$u = \begin{pmatrix} v \\ \omega \end{pmatrix}$

The input vector u is constrained as follows: ∥v∥≦V_(max) ∥ω∥≦ω_(max) where V_(max) and ω_(max) are maxima of the vehicle's speed and angular velocity respectively.

In this embodiment, the state vector x of the vehicle 2 is linearised by replacing the heading 4 with components of the speed v in the x- and y-directions, i.e.:

$x = \begin{pmatrix} r_{x} \\ r_{y} \\ {v\;\sin\;\theta} \\ {v\;\cos\;\theta} \end{pmatrix}$

The dynamics of the system are given by a linear time-invariant system of the form: {dot over (x)}=Ax+Bu

Because the vehicle's speed v is part of the state vector x, the linearised form of the input vector u comprises acceleration components in the x- and y-directions, i.e.

$u = \begin{pmatrix} a_{x} \\ a_{y} \end{pmatrix}$

FIG. 2 is a schematic illustration (not to scale) of a scenario in which the vehicle 2 travels along a trajectory 6 (indicated by a dotted line in FIG. 2).

The trajectory 6 is divided into a series of points p₀, p₁, p₂, . . . , p_(N). The ith point on the trajectory 6, i.e. p_(i), is a position occupied by the vehicle at an ith time-step.

State vectors of the vehicle 2 at each of the points p₀, . . . , p_(N) of the trajectory 6 are x(0), . . . , x(N) respectively. In other words, x(i) denotes the state of the vehicle 2 at the ith time-step. Also, an initial state of the vehicle 2, i.e. a state of the vehicle 2 at the start of the trajectory 6, is x₀.

A desired terminal position of the vehicle 2 is indicated in FIG. 2 by the point P.

In this embodiment, the starting position p₀ and the desired terminal position P are identified for the vehicle 2, for example by a user/operator of the vehicle 2, by any appropriate manner. As used herein, the terminology “identified” includes any appropriate form of identifying, selecting, choosing, establishing, acquiring etc.

A state of the vehicle 2 at the desired terminal position P is X_(P).

Thus, it is desirable for the vehicle 2 to follow a trajectory 6 such that the distance d between the desired terminal position P, and the terminal position of the vehicle 2 after following the trajectory 6, i.e. the point p_(N), is minimised. In other embodiments, a ‘stage cost’ is also minimised.

In this embodiment, such an optimal trajectory (a trajectory that minimises the distance between P and p_(N)) is determined by a trajectory planner (not shown in the Figures).

In this embodiment, the determination of the optimal trajectory is implemented using the following equation:

$\min\limits_{u,x}{g\left( {x,u} \right)}$ where the value function is given by g(x,u)=∥[I 0](x(N)−x _(P))∥ In other embodiments this may also include cost term associated with traversing the trajectory, such as time or distance.

In this embodiment the optimisation problem is subject to the following constraints: x(0)=x ₀ x(k+1)=Ax(k)+Bu(k) ([0 I]x(k),u(k))εL where: g(x,u) is a value function that defines optimality;

-   -   I is a (2×2) identity matrix; and     -   L is a net of acceptable states (v_(x),v_(y)) and controls         (a_(x), a_(y)). Here, v_(x) is the speed of the vehicle 2 in the         x-direction, v_(y) is the speed of the vehicle 2 in the         y-direction, a_(x) is the acceleration of the vehicle 2 in the         x-direction, and a_(y) is the acceleration of the vehicle 2 in         the y-direction.

In this embodiment, L is used to constrain the magnitude of the speed and the acceleration of the vehicle 2.

Thus, the above equation provides that the Euclidean distance between x(N) and xp is minimised (the velocity components of the state vector, x, are multiplied by 0 to remove them).

In this embodiment, the above described constraint equations are used. However, in other embodiments different constraint equations may be used instead of, or in addition to, some or all of the above constraint equations. For example, in other embodiments, constraint equations may be used that provide that a control effort for the vehicle 2, or a number of time-steps to reach the destination, is minimised.

The above equation for determining the optimal trajectory does not take into account the heading 4 of the vehicle 2. Thus, a trajectory determined using this equation alone may, at a certain point, include motion that cannot be performed by the vehicle 2 used in this embodiment. For example, the trajectory 6 determined as described above may require that, at a certain point, the vehicle 2 travels in a direction that is perpendicular to the vehicle's heading 4 at that point. Such a trajectory cannot be followed by the vehicle 2 of this embodiment (i.e. a land-based vehicle) because it would require an infinitely large turn-rate.

Conventionally, this constraint on the vehicle's turn rate is treated as an acceleration constraint. “Receding Horizon Control In Unknown Environments: Experimental results”, Markus Deittert, Arthur Richards, and George Mathews, ICRA, Achorage, Ak., USA, May 2010, which is incorporated herein by reference, shows an implementation in which the magnitude of the vehicle's input vector u (i.e. an acceleration) is limited in relation to the vehicles maximum velocity, V_(max), such that a minimum turn radius, R_(min), is enforced, i.e.:

$a_{{ma}\; x} = \frac{V_{m\;{ax}}^{2}}{R_{m\; i\; n}}$

Thus, conventionally, trajectories that include a turn having a turn radius of less than R_(min) tend to be avoided by the vehicle travelling at velocities close to its maximum, V_(max).

In this embodiment, an acceleration of the vehicle 2 that changes the direction of the vehicle 2, but that does not significantly change the norm of the vehicle's velocity is used.

This is achieved by requiring that the acceleration primarily acts in directions that are normal (i.e. perpendicular) to the vehicle's velocity vector v. In this embodiment, the vehicle's velocity vector at the kth time step, is: v(k)=(v _(x) ,v _(y))^(T) where: v_(x) is the speed of the vehicle 2 in the x-direction at the kth time-step; and

-   -   v_(y) is the speed of the vehicle 2 in the y-direction at the         kth time-step.

FIG. 3 is a schematic illustration that is useful for understanding an optional additional technique that may be implemented with this embodiment, and will be described in more detail later below.

FIG. 4 is a process flow chart showing certain steps of a process of determining a constrained acceleration vector for the vehicle 2 according to the first additional technique.

FIG. 5 is a schematic illustration that is useful in the understanding of the process of FIG. 4.

At step s2, the velocity vector v of the vehicle 2 is normalised.

At step s4, the normalised velocity vector {circumflex over (v)} is multiplied by a scalar quantity λ.

At step s6, two convex polynomial approximations of a circle, hereinafter referred to as “the first polygon” and the “second polygon” and indicated in FIG. 5 by the reference numerals 14 and 16 respectively, are determined.

The first polygon 14 is centred at a point λ{circumflex over (v)}.

The second polygon 16 is centred at a point −λ{circumflex over (v)}.

The radii of the first and second polygons 14, 16 are indicated in FIG. 5 by dotted arrows. These radii are larger than the length of the vector λ{circumflex over (v)}. This provides that the first and second polygons 14, 16 overlap to some extent.

At step s8, an acceleration vector a is determined such that it lies within an overlap 18 of the first polygon 14 with the second polygon. This provides that the acceleration applied to the vehicle 2 may be relatively large when acting in a direction that is substantially perpendicular to the velocity vector v of the vehicle 2, but is relatively small when acting in a direction that is substantially parallel to the velocity vector v of the vehicle 2.

If λ is selected to be a relatively large value, e.g. λ=10, the resulting overlap 18 is relatively small in a direction that is substantially parallel to the velocity of the vehicle v, but is relatively large in a direction that is substantially normal to the velocity vector v. The value of λ may advantageously be selected depending on the application.

Due to the first and second polygons 14, 16 being centred around points, the position of which depends on the normalised velocity vector {circumflex over (v)}, if the vehicle 2 slows down, the overlap 18 increases relative to the velocity vector {circumflex over (v)}. This advantageously tends to provide that, at low speeds, a relatively large acceleration may be applied to the vehicle in a direction that is substantially parallel to the velocity vector v, thereby allowing the magnitude of the velocity vector v to be increased.

The process of determining an acceleration vector for the vehicle 2 described above with reference to FIGS. 4 and 5 is equivalent to determining an acceleration vector

${a(k)} = \begin{pmatrix} a_{x} \\ a_{y} \end{pmatrix}$ for the vehicle 2 that is constrained as follows:

  ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circs)] ${{\left( {\frac{{a_{x}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{x}(k)}}{V_{{ma}\; x}}} \right){\cos\left( \frac{2\pi\; m}{N_{{circs}\;}} \right)}} + {\left( {\frac{{a_{y}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{y}(k)}}{V_{{ma}\; x}}} \right){\sin\left( \frac{2\pi\; m}{N_{circs}} \right)}}} \leq \sqrt{1 + \lambda^{2}}$ where: v(k)=(v _(x) ,v _(y))^(T); and

N_(circs) is the number of corners of the first and second polygons 14, 16. In other embodiments, the first polygon 14 may have a different number of corners to the second polygon 16 (and the above equation is modified accordingly).

Thus, a technique by which a vehicle trajectory may be determined using a trajectory planner is provided. The above described constraints applied to the trajectory planner, i.e. the constraints on the vector v, advantageously tend to provide that the determined trajectory is able to be followed by a vehicle (e.g. a wheeled land vehicle) that has a curvature limit.

The above described embodiment, in which a constraint on the acceleration vector a of the vehicle 2 is applied, advantageously tends to provide that the trajectory planner tends not to be able to produce trajectories which require, at a particular point in time, the vehicle 2 to travel perpendicular to its heading 4.

A further advantage provided by the above described trajectory planner constraint is that the performance of the trajectory planner and/or the vehicle, in particular when the vehicle travels at relatively low speeds, tends to be improved compared to conventional approaches.

A further advantage provided by the above described embodiment is that a linear approximation of the vehicle's dynamics is advantageously constrained. In particular, in the above embodiment the magnitude of the vehicle acceleration in a direction perpendicular to the vehicle's heading is constrained. This is achieved by scaling the components of the acceleration vector (i.e. the components a_(x) and a_(y)) depending of the respective components of the velocity vector (i.e. the components v_(x) and v_(y) respectively). In other words, constraints on the trajectory of the vehicle are implemented in the acceleration space of the vehicle. In particular, the magnitude and the direction of the vehicle's acceleration is constrained.

In other embodiments, the above described approach may be combined with one or more of the following additional optional techniques.

A first additional technique involves constraining a trajectory planner such that it may not plan trajectories that require the vehicle 2 to travel in regions that limitations on the vehicle's turning circle prevent it from travelling in.

In particular, for curvature limited vehicles, which cannot turn on the spot, there exists a circular area to each side of the vehicle that cannot be reached by turning directly into it, for example without repeatedly reversing and advancing. The radius of these circles is substantially equal to the vehicle's minimum turn radius, R_(min).

FIG. 6 is a schematic illustration showing the vehicle 2 and two regions, hereinafter referred to as the “first region 20” and the “second region 22”, into which the vehicle 2 may not directly move, i.e. turn directly into.

The vehicle 2 may not move into the first or second region due to the vehicle's curvature limit. The first and second region 20, 22 are avoided by the trajectory planner when planning a trajectory. The first and second regions 20, 22 may be considered to be ‘obstacles’ that are to be avoided when planning a trajectory of the vehicle 2.

In embodiments in which the first additional technique is implemented, the first and second regions 20, 22 are approximated by polygons, and the constraint on the trajectory of the vehicle 2 is implemented as follows.

A unit vector normal to the initial speed vector of the vehicle v(0) is:

${\hat{v}}_{n} = {\left( {v_{nx},v_{ny}} \right)^{T} = \frac{\left( {{v_{y}(0)},{- {v_{x}(0)}}} \right)}{\sqrt{{v_{x}^{2}(0)} + {v_{y}^{2}(0)}}}}$ where: v_(nx) the component of the vector {circumflex over (v)}_(n) in the x-direction; and

-   -   v_(ny) is the component of the vector {circumflex over (v)}_(n)         in the y-direction.

An initial position of the vehicle 2 is given by: r₀=(r_(0x),r_(0y))^(T).

In embodiments in which the first additional technique is implemented, the constraint on the trajectory of the vehicle 2 supplied by one of the region 20, 22 is implemented as follows:

  ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(regions)] ${{\left( {{r_{x}(k)} - r_{0x} - {R_{m\; i\; n}v_{nx}}} \right){\cos\left( \frac{2\pi\; m}{N_{regions}} \right)}} + {\left( {{r_{y}(k)} - r_{0y} - {R_{m\; i\; n}v_{ny}}} \right){\sin\left( \frac{2\pi\; n}{N_{regions}} \right)}}} \geq {R_{m\; i\; n} - {C_{L}{p\left( {k,m} \right)}}}$ with

∀k ∈ [0, …  , N_(steps)] ${\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{regions} - 1}$ where: N_(steps) is the number of time-steps;

-   -   N_(regions) is the number of corners of the polygons used to         approximate the first or second regions 20, 22. In this         embodiment, the polygons used to approximate the first and         second regions 20, 22 comprise the same number of points.         However, in other embodiments these polygons may comprise a         different number of points. In such cases the above equations         may be modified accordingly;     -   p(k,m) is a binary decision variable. In this embodiment p(k,m)         is a matrix of size N_(steps)×N_(regions) with an entry (equal         to zero or one) for each time-step k and polygon corner m. In         this embodiment, p(k,m) is used to deactivate constraints; and

C_(L) is an arbitrary constant greater than R_(min). In this embodiment C_(L) is equal to 2×R_(min).

The constant C_(L) may be advantageously selected depending on the application.

Constraining the trajectory in this way advantageously tends to provide that vehicle 2 may follow the trajectory, even at relatively low speeds. This tends to be particularly useful when planning a trajectory from a resting position or in a cluttered surrounding.

Thus, a first additional, optional technique by which a vehicle trajectory may be determined using a trajectory planner is provided. The above described constraints applied to the trajectory planner, i.e. the constraints on a position vector r(k)=(r_(x),r_(y)) of the vehicle, advantageously tend to provide that the determined trajectory is able to be followed by a vehicle (e.g. a wheeled land vehicle) that has a curvature limit.

The above described embodiment, in which the trajectory is required to provide that the vehicle 2 avoids regions close to, and either side of, the vehicle (as described above with reference to FIG. 6), advantageously tends to provide that the trajectory planner tends not to be able to produce trajectories which require, at a particular point in time, the vehicle 2 to travel perpendicular to its heading 4.

A further advantage provided by the above described trajectory planner constraint is that the performance of the trajectory planner and/or the vehicle, in particular when the vehicle travels at relatively low speeds, tends to be improved.

A further advantage provided by the above described embodiment is that a linear approximation of the vehicle's dynamics is advantageously constrained. In particular, in the above embodiment the magnitude of the vehicle acceleration in a direction perpendicular to the vehicle's heading is constrained. This is achieved by determining the components of the position vector (i.e. the components r_(x) and r_(y)), which define a point that the vehicle will be moved to at a particular point in time (in effect the heading of the vehicle at a point in time) depending of the respective components of the velocity vector (i.e. the components v_(x) and v_(y) respectively). In other words, constraints on the trajectory of the vehicle are implemented in the positional space (i.e. x-y space) of the vehicle.

The above described conventional approach for transforming the vehicle's turn rate limit into an acceleration constraint tends to fail at low speeds. In particular, if the vehicle slows down, the enforcement of

$a_{{ma}\; x} = \frac{V_{{ma}\; x}^{2}}{R_{m\; i\; n}}$ tends to result in trajectories comprising turns with a turn radius less than R_(min).

A second additional technique involves imposing a minimum constraint to the vehicle's velocity. This constraint provides that the minimum speed of the vehicle 2, denoted hereinafter as “V_(min)”, is close to the maximum speed of the vehicle 2.

In this embodiment, this minimum speed constraint is enforced as follows.

FIG. 3 is a schematic illustration representing a space of the vehicle's velocity, hereinafter referred to as the “velocity plane 8”.

A convex polynomial approximation of a circle is indicated by the reference numeral 10 and is hereinafter referred to as the “polygon”. The polygon 10 approximates a circle having a radius equal to V_(min), indicated in FIG. 3 by an arrow and the reference numeral 12. Also, the polygon 10 is centred at the origin ◯′ of the velocity plane 8.

The velocity vector v(k) is constrained to remain outside the polygon 10, i.e.

∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circ)] ${{{v_{x}(k)}{\sin\left( \frac{2\pi\; m}{N_{circ}} \right)}} + {{v_{y}(k)}{\cos\left( \frac{2\pi\; m}{N_{circ}} \right)}}} \geq {V_{m\; i\; n} - {{p\left( {k,m} \right)} \cdot C_{L}}}$ with ∀k ∈ [0, …  , N_(steps)] ${\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{circ} - 1}$ Where: N_(steps) is the number of time-steps;

-   -   N_(circ) is the number of corners of the polygon 10;     -   p(k,m) is a binary decision variable. In this embodiment p(k,m)         is a matrix of size N_(steps)×N_(circ) with an entry (equal to         zero or one) for each time-step k and polygon corner m. In this         embodiment, p(k,m) is used to deactivate constraints; and     -   C_(L) is an arbitrary constant. In this embodiment C_(L) is         relatively large compared to V_(min).

The above velocity constraint is satisfied either by v(k) being outside the polygon 10, or by each value in the matrix p(k,m) being equal to one. The constraint

${\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{circ} - 1}$ provides that only N_(circ)−1 entries of p(k,m) are equal to one. Thus, v(k) must be outside of the polygon 10.

Thus, a second additional, optional technique by which a vehicle trajectory may be determined using a trajectory planner is provided. The above described constraints applied to the trajectory planner, i.e. the constraints on the vector v, advantageously tend to provide that the determined trajectory is able to be followed by a vehicle (e.g. a wheeled land vehicle) that has a curvature limit.

The above described embodiment, in which a constraint on the velocity vector v of the vehicle 2 is applied, advantageously tends to provide that the trajectory planner tends not to be able to produce trajectories which require, at a particular point in time, the vehicle 2 to travel perpendicular to its heading 4.

A further advantage provided by the second additional optional trajectory planner constraint is that the performance of the trajectory planner and/or the vehicle, in particular when the vehicle travels at relatively low speeds, tends to be improved compared to the conventional approach. Performance may, for example, be measured as the error between the turn radii of the optimal linear trajectory and the turn radius limit of the non-linear vehicle dynamics.

A further advantage provided by the above described embodiment is that a linear approximation of the vehicle's dynamics is advantageously constrained. In particular, in the above embodiment the magnitude of the vehicle acceleration in a direction perpendicular to the vehicle's heading is constrained. This is achieved by constraining the velocity vector of the vehicle as described above. In other words, constraints on the trajectory of the vehicle are implemented in the velocity space of the vehicle.

A third additional technique involves linearly approximating the maximum acceleration of the vehicle 2.

The magnitude of the velocity of the vehicle |v|, and maximum acceleration of the vehicle a_(max), are related by a nonlinear function. Conventionally, this non-linear function cannot be used directly in the optimisation of the trajectory 6. However, in this embodiment the nonlinear function is approximated by a collection of linear functions.

In embodiments in which the third additional technique is implemented the ∞-norm of the velocity of the vehicle v(k)=(v_(x),v_(y))^(T) is used. This advantageously relates the amount of acceleration available to the trajectory planner to the largest component within the speed vector, v.

FIG. 7 is a schematic illustration (not to scale) of a graph showing the non-linear relationship between maximum acceleration of the vehicle a_(max), the ∞-norm of the velocity of the vehicle |v|, and the minimum turn radius R_(min). In FIG. 7 this non-linear function is indicated by the reference numeral 24, and each of the collection of linear functions that are used to approximate the nonlinear function 24 are indicated by the reference numeral 26.

In embodiments in which the third additional technique is implemented, the following constraints on the acceleration vector of the vehicle a(k)=(a_(x),a_(y))^(T) and the velocity vector of the vehicle v(k)=(v_(x),v_(y))^(T) are implemented:

$\begin{matrix} {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{m \in \left\lbrack {1,\ldots\mspace{14mu},N_{cir}} \right\rbrack}}}{{{{a_{x}(k)}{\sin\left( \frac{2\pi\; m}{N_{cir}} \right)}} + {{a_{y}(k)}{\cos\left( \frac{2\pi\; m}{N_{cir}} \right)}}} \leq {a_{{ma}\; x}(k)}}} & (i) \\ {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{p \in \left\lbrack {1,\ldots\mspace{14mu},N_{pol}} \right\rbrack}}}{{{{m(p)}{v_{x}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 1\left( {k,p} \right)}}}{{{{m(p)}{v_{y}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 2\left( {k,p} \right)}}}} & ({ii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{\sum\limits_{p = 1}^{N_{pol}}\left( {{b\; 1\left( {k,p} \right)} + {b\; 2\left( {k,p} \right)}} \right)} \leq {{2N_{pol}} - 1}}} & ({iii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{a_{{ma}\; x}(k)} \geq 0}} & ({iv}) \end{matrix}$ where: N_(steps) is a number of time-steps;

-   -   N_(cir) is a number of corners of a polygon (a convex polynomial         approximation of a circle). This polygon approximates a circle         having a radius equal to a_(max);     -   N_(pol) is the number of linear functions 26 that are used to         approximate the nonlinear function 24;     -   m(p) is the value of the gradient of the pth linear function 26     -   c(p) is the value of the ‘a_(max)-intercept’ of the pth linear         function 26, i.e. the value of the pth linear function 26 if the         value of the velocity was zero;     -   b1(k,p) and b2(k,p) are constants the values of which provide         that at least one of the equations of the second constraint         (i.e. constraint (ii)) are fulfilled by a value of a_(max); and     -   C_(L) is a constant.

The first of the above constraints, i.e. constraint (i), provides that the magnitude of the acceleration vector a(k)=(a_(x),a_(y))^(T) is limited, i.e. that the magnitude of the acceleration vector is less than or equal to a_(max) for a particular time-step.

The second of the above constraints, i.e. constraint (ii), provides that the maximum acceleration of the vehicle a_(max) is less than or equal to the approximated non-linear function 24;

The third of the above constraints, i.e. constraint (iii), provides that a trajectory planner implementing the above constraints (i)-(iv) bases a value of a_(max) on the x-component or y-component of the velocity vector v. The ∞-norm of the velocity vector v is thereby implemented.

The fourth of the above constraints, i.e. constraint (iv), provides that the maximum acceleration a_(max) is always positive.

Thus, a third additional, optional technique by which a vehicle trajectory may be determined using a trajectory planner is provided. The above described constraints applied to the trajectory planner, i.e. that the magnitude of acceleration of the vehicle (not the direction) is constrained, advantageously tend to provide that the determined trajectory is able to be followed by a vehicle (e.g. a wheeled land vehicle) that has a curvature limit.

The above described embodiment, in which the trajectory is constrained in such a way that inter alia the maximum value of the acceleration of the vehicle a_(max) is a function of either the x-component or y-component of the velocity vector v, advantageously tends to provide that the trajectory planner tends not to be able to produce trajectories which require, at a particular point in time, the vehicle 2 to travel perpendicular to its heading 4.

A further advantage provided by the above described trajectory planner constraint is that the performance of the trajectory planner and/or the vehicle, in particular when the vehicle travels at relatively low speeds, tends to be improved.

A further advantage provided by the above described embodiment is that a linear approximation of the vehicle's dynamics is advantageously constrained. In particular, in the above embodiment the magnitude of the vehicle acceleration in a direction perpendicular to the vehicle's heading is constrained. This is achieved by linear approximation of the acceleration limit of the vehicle (i.e. the maximum acceleration). In other words, constraints on the trajectory of the vehicle are implemented in the acceleration space of the vehicle. In particular, the magnitude of the maximum acceleration is constrained depending on the velocity vector v.

One or more of the above described optional additional techniques for constraining a trajectory determined by a trajectory planner may advantageously be implemented in conjunction with, or instead of, the acceleration constraints described above with reference to FIGS. 4 and 5. It may be desirable to use different sets of constraints in different situations.

Different sets of constraints (i.e. the acceleration constraints described above with respect to FIGS. 4 and 5, or the constraints relating to any of the above described additional technique) may be implemented by a trajectory planner contemporaneously or in series. How the constraints are implemented, i.e. in what order and for what proportion of the vehicle's trajectory, may advantageously be selected depending on the application. For example, in another embodiment two sets of constraints may be applied at different times, e.g. the positional constraints of the first additional technique may be implemented for a first half of a vehicle's trajectory, and a different constraint may be implemented for the later half. In another embodiment, a point in time at which a constraint set being implemented by a trajectory planner is changed may be decided using a binary decision variable. In such a way, the change over from one set of constraints to another can advantageously be tied to the vehicle's state.

Apparatus, including the trajectory planner (not shown in the Figures), for implementing the above arrangement, and performing any of the above described method steps, may be provided by configuring or adapting any suitable apparatus, for example one or more computers or other processing apparatus or processors, and/or providing additional modules. The apparatus may comprise a computer, a network of computers, or one or more processors, for implementing instructions and using data, including instructions and data in the form of a computer program or plurality of computer programs stored in or on a machine readable storage medium such as computer memory, a computer disk, ROM, PROM etc., or any combination of these or other storage media.

It should be noted that certain of the process steps depicted in the flowchart of FIG. 4 and described above may be omitted or such process steps may be performed in differing order to that presented above and shown in FIG. 4. Furthermore, although all the process steps have, for convenience and ease of understanding, been depicted as discrete temporally-sequential steps, nevertheless some of the process steps may in fact be performed simultaneously or at least overlapping to some extent temporally.

In the above embodiments, the vehicle is a land-based vehicle, e.g. a vehicle comprises wheels and/or tracks. However, in other embodiments the vehicle may be any appropriate vehicle that has a curvature limit, e.g. a boat, submarine, or amphibious vehicle. Also, the vehicle may be manned or unmanned.

In the above embodiments, the constraints that are implemented by the trajectory planner are expressed by the relevant above described equations. In other embodiments, one or more of the constraints may be implemented using a different appropriate equation so as to provide an equivalent constraining effect on the trajectory planner, and/or provide the equivalent functionality to that described above.

In the above embodiments, the polygons used in the linear approximations may comprise any appropriate number of corners. Generally, the greater the number of corners used for the polygon(s), the greater the accuracy of the approximation to the circles it/they represent tend to be. However, the greater the number of corners of the polygon(s), the more processing power is required. Thus, a trade-off exists between accuracy and processing power. The number of corners for each of the polygons used may be selected. A trade-off exists between the accuracy of the approximation and the processing power required to perform the approximation in a certain amount of time. The number of corners for each polygon may be advantageously selected to achieve a desired balance between accuracy and processing power. 

The invention claimed is:
 1. A method for controlling a land vehicle having a curvature limit according to a trajectory, the method comprising: identifying a starting position (p₀) for the land vehicle; identifying a desired terminal position (P) for the land vehicle; defining a linear approximation of by one or more processors, dynamics of the land vehicle with a state vector (x) that indicates land vehicle velocity with linearized variables that lack non-linear heading information; determining, based on the starting position (p₀), the desired terminal position (P), and the linear approximation, by the one or more processors, the trajectory for the land vehicle that optimizes a linear time invariant system having a form {dot over (x)}=Ax+Bu against a cost function and with a constraint applied to the state vector (x) to determine an acceleration input (u) to be applied by the land vehicle to accomplish the trajectory; and controlling the land vehicle to follow the trajectory; wherein to prevent the trajectory from including curvatures below the curvature limit of the land vehicle, the constraint requires that the acceleration is larger when acting in a direction that is substantially perpendicular to the velocity of the land vehicle than when acting in a direction that is substantially parallel to the velocity of the vehicle is implemented by calculating:   ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circs)] ${{\left( {\frac{{a_{x}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{x}(k)}}{V_{{ma}\; x}}} \right){\cos\left( \frac{2\pi\; m}{N_{circs}} \right)}} + {\left( {\frac{{a_{y}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{y}(k)}}{V_{{ma}\; x}}} \right){\sin\left( \frac{2\pi\; m}{N_{circs}} \right)}}} \leq \sqrt{1 + \lambda^{2}}$ where: v_(x) (k) is the velocity of the land vehicle in an x-direction; v_(y) (k) is the velocity of the land vehicle in a y-direction; V_(max) is a maximum velocity of the land vehicle; N_(steps) is a number of time-steps for the trajectory; N_(circs) is a number of corners of a polygon; R_(min) is a minimum turn radius of the land vehicle; and λ is a constant wherein the linear approximation of the vehicle dynamics is further constrained by specifications that a magnitude of a velocity of the vehicle is greater than or equal to a threshold value for the velocity.
 2. A method according to claim 1, wherein the constraint that a magnitude of the velocity of the vehicle is greater than or equal to a threshold value for the velocity is implemented by calculating: and ∀k ∈ [0, …  , N_(steps)] ${\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{circ} - 1}$ where: v_(x)(k) is the velocity of the vehicle in an x-direction; v_(y)(k) is the velocity of the vehicle in a y-direction; V_(min) is the threshold velocity value; N_(steps) is a number of time-steps for the trajectory; N_(circ) is a number of corners of a polygon; p(k,m) is a binary decision variable; and C_(l) is a constant.
 3. A method according to claim 2, wherein C_(l) is relatively large compared to V_(min).
 4. A method according to claim 1, wherein: the linear approximation is further constrained by a specification that the vehicle not travel in a given region; a centre point of the region is at a distance from a predetermined point on the vehicle substantially equal to a minimum turn radius of the vehicle, in a direction substantially perpendicular to a velocity of the vehicle; and a distance from the centre point of the region to a point on a perimeter of the region is greater than or substantially equal to the minimum turn radius of the vehicle.
 5. A method according to claim 4, wherein the constraint that the vehicle not travel in a given region is implemented by calculating:   ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(regions)] ${{\left( {{r_{x}(k)} - r_{0x} - {R_{m\; i\; n}v_{nx}}} \right){\cos\left( \frac{2\pi\; m}{N_{regions}} \right)}} + {\left( {{r_{y}(k)} - r_{0y} - {R_{m\; i\; n}v_{ny}}} \right){\sin\left( \frac{2\pi\; m}{N_{regions}} \right)}}} \geq {R_{m\; i\; n} - {C_{2}{p\left( {k,m} \right)}}}$   with   ∀k ∈ [0, …  , N_(steps)] $\mspace{20mu}{{\sum\limits_{m = 1}^{N_{circ}}{p\left( {k,m} \right)}} = {N_{regions} - 1}}$ where: v_(nx) is a normalised velocity of the vehicle in an x-direction; v_(ny) is a normalised velocity of the vehicle in a y-direction; r_(x)(k) is a position of the vehicle from an origin in the x-direction; r_(y)(k) is a position of the vehicle from an origin in the y-direction; r_(0x) is an initial position of the vehicle from an origin in the x-direction; r_(0y) is an initial position of the vehicle from an origin in the y-'direction; N_(steps) is a number of time-steps for the trajectory; N_(regions) is a number of corners of a polygon; R_(min) is a minimum turn radius of the vehicle; p(k,m) is a binary decision variable; and C₂ is a constant.
 6. A method according to claim 1, wherein the linear approximation is further constrained by specifications that: an acceleration of the vehicle during the trajectory is less than a threshold value for the acceleration; and the threshold value for the acceleration is dependent on an infinity norm of a velocity of the vehicle.
 7. A method according to claim 6, wherein the constraint that the linear approximation is further constrained by specifications that an acceleration of the vehicle during the trajectory is less than a threshold value for the acceleration, and the threshold value for the acceleration is dependent on an infinity norm of a velocity of the vehicle is implemented by calculating: $\begin{matrix} {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{m \in \left\lbrack {1,\ldots\mspace{14mu},N_{cir}} \right\rbrack}}}{{{{a_{x}(k)}{\sin\left( \frac{2\pi\; m}{N_{cir}} \right)}} + {{a_{y}(k)}{\cos\left( \frac{2\pi\; m}{N_{cir}} \right)}}} \leq {a_{m\;{ax}}(k)}}} & (i) \\ {{{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}},{\forall{p \in \left\lbrack {1,\ldots\mspace{14mu},N_{pol}} \right\rbrack}}}{{{{m(p)}{v_{x}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 1\left( {k,p} \right)}}}{{{{m(p)}{v_{y}(k)}} + {c(p)}} \geq {{a_{{ma}\; x}(p)} - {C_{L}b\; 2\left( {k,p} \right)}}}} & ({ii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{\sum\limits_{p = 1}^{N_{pol}}\left( {{b\; 1\left( {k,p} \right)} + {b\; 2\left( {k,p} \right)}} \right)} \leq {{2N_{pol}} - 1}}} & ({iii}) \\ {{\forall{k \in \left\lbrack {0,\ldots\mspace{14mu},N_{steps}} \right\rbrack}}{{a_{{ma}\; x}(k)} \geq 0}} & ({iv}) \end{matrix}$ where: a_(max) is a maximum acceleration of the vehicle; a_(x) is an acceleration of the vehicle in an x-direction; a_(y) is an acceleration of the vehicle in a y-direction; v_(x) is a velocity of the vehicle in an x-direction; v_(y) is a velocity of the vehicle in a y-direction; N^(steps) is a number of time-steps for the trajectory; N_(cir) is a number of corners of a polygon; N_(pol) is a number of linear functions, each linear function relating the maximum acceleration of the vehicle to an infinity norm of the velocity of the vehicle; m(p) is a value of a gradient of a pth linear function; c(p) is a value of pth linear function if the value of the velocity of the vehicle is zero; b1(k, p) and b2(k, p) are constants; and C_(L) is a constant.
 8. A method according to claim 1, comprising: determining, by the one or more processors, one or more further linear approximations of the dynamics of the vehicle; and determining, based on the one or more further linear approximations, by the one or more processors, one or more further trajectories for the vehicle; wherein the trajectory and the one or more further trajectories are for implementation by the vehicle in series; each of the one or more further linear approximations is constrained by one or more of the following: (i) a specification that a magnitude of a velocity of the vehicle is greater than or equal to a threshold value for the velocity; (ii) specifications that: an acceleration applied to the vehicle at a point on the trajectory has a first value when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; and an acceleration applied to the vehicle at a point on the trajectory is smaller than the first value when the acceleration acts in a direction that is substantially perpendicular to the velocity of the vehicle; (iii) a specification that the vehicle not travel in a given region; a centre point of the region is at a distance from a predetermined point on the vehicle substantially equal to a minimum turn radius of the vehicle, in a direction substantially perpendicular to a velocity of the vehicle; and a distance from the centre point of the region to a point on a perimeter of the region is greater than or substantially equal to the minimum turn radius of the vehicle; and (iv) specifications that an acceleration of the vehicle during the trajectory is less than a threshold value for the acceleration, and the threshold value for the acceleration is dependent on an infinity norm of a velocity of the vehicle.
 9. Apparatus comprising one or more processors and at least one non-transitory computer readable storage medium having encoded thereon executable instructions that when executed by the one or more processors, cause the one or more processors to carry out a method for controlling a land vehicle having. a curvature limit, the method comprising: identifying a starting position (p₀) for the land vehicle; identifying a desire terminal position (P) for the land vehicle; defining a linear approximation of by one or more processors, dynamics of the land vehicle with a state vector (x) that indicates and vehicle velocity with linearized variables that lack non-linear headind information; determining, based on the starting position (p₀), the desired terminal position (P), and the linear approximation, by the one or more processors, the trajectory for the land vehicle that optimizes a linear time invariant system having a form {dot over (x)}=Ax+Bu against a cost function and with a constraint applied to the state vector (x) to determine an acceleration input (u) to be applied by the land vehicle to accomplish the trajectory, and controlling the land vehicle to follow the trajectory; wherein to prevent the trajectory from including curvatures below the curvature limit of the land vehicle, the constraint requires that the acceleration is larger when acting in a direction that is substantially perpendicular to the velocity of the land vehicle than when acting in a direction that is substantially parallel to the velocity of the vehicle is implemented by calculating:   ∀k ∈ [0, …  , N_(steps)], ∀m ∈ [1, …  , N_(circs)] ${{\left( {\frac{{a_{x}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{x}(k)}}{V_{{ma}\; x}}} \right){\cos\left( \frac{2\pi\; m}{N_{circs}} \right)}} + {\left( {\frac{{a_{y}(k)}R_{m\; i\; n}}{V_{{ma}\; x}} - \frac{\lambda\;{v_{y}(k)}}{V_{{ma}\; x}}} \right){\sin\left( \frac{2\pi\; m}{N_{circs}} \right)}}} \leq \sqrt{1 + \lambda^{2}}$ where: v_(x) (k) is the velocity of the land vehicle in an x-direction; v_(y)(k) is the velocity of the land vehicle in a v-direction; V_(max) is a maximum velocity of the land vehicle; N_(steps) is a number of time-steps for the trajectory; N_(circs) is a number of corners of a polygon; R_(min) is a minimum turn radius of the land vehicle: and λ is a constant wherein the linear approximation of the vehicle dynamics is further constrained by specifications that a magnitude of a velocity of the vehicle is greater than or equal to a threshold value for the velocity. 